Nonlinear Whistlerons* 



^1" 
O 

o 

(N 
•*-> 
O 

O 

(N 



43 



B. Eliasson, I. Kourakis and P. K. Shukla 

Institut fur Theoretische Physik IV, Fakultat fur Physik und Astronomie, 
Ruhr-Universitat Bochum, D-44780 Bochum, Germany 
(Dated: Received 12 October 2004) 

Recently, observations from laboratory experiments have revealed amplitude modulation of 
whistlers by low-frequency perturbations. We here present theoretical and simulation studies of 
amplitude modulated whistler solitary waves (whistlerons) and their interaction with background 
low-frequency density perturbations created by the whistler ponderomotive force. We derive a non- 
linear a nonlinear Schrodinger equation which governs the evolution of whistlers in the presence of 
finite-amplitude density perturbations, and a set of equations for arbitrary large amplitude density 
perturbations in the presence of the whistler ponderomotive force. The governing equations stud- 
ied analytically in the small amplitude limit, and are solved numerically to show the existence of 
large scale density perturbations that are self-consistently created by localized whistlerons. Our 
numerical results are in good agreement with recent experimental results where the the formation 
of modulated whistlers and solitary whister waves were formed. 
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II. DERIVATION OF THE GOVERNING 
EQUATIONS 
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Almost three decades ago, Stenzel pj experimentally 
demonstrated the creation of a magnetic field-aligned 
density cavities by the ponderomotive force of localized 
electron whistlers. Observations from a recent labora- 
tory experiment exhibit the creation of modulated 
whistler wavepackets due to nonlinear effects. Further- 
more, instruments on board the CLUSTER spacecraft 
have been observing broadband intense electromagnetic 
waves, correlated density fluctuations and solitary waves 
near the Earth's plasmapause, magnetopause and fore- 
shock revealing signatures of whistler turbulence in 
the presence of density depletions and enhancements. 
The Freja satellite |j] also observed the formation of enve- 
lope whistler solitary waves correlated with density cav- 
ities in the plasma. 

A theoretical investigation has in the past predicted 
the self- channeling of electron whistlers and the creation 
of a localized density hump ,5] . Taking into account the 
spatio-temporal dependent whistler ponderomotive force 
[a, Q, investigations of the modulation and filamenta- 
tion of finite amplitude whistlers interacting with mag- 
netosonic waves [g, El an d ion-acoustic perturbations 
have been carried out. 

In this article, we investigate nonlinearly interacting 
electron whistlers and arbitrary large amplitude ion- 
acoustic perturbations, by using computer simulations, 
and we find analytical expressions for whistlerons in the 
low-amplitude limit 
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Let us consider the propagation of nonlinearly coupled 
whistlers and ion-acoustic perturbations in a fully ionized 
electron-ion plasma in a uniform external magnetic field 
zBq , where z is the unit vector along the z direction and 
Bo is the magnitude of the magnetic field strength. We 
consider the propagation of right-hand circularly polar- 
ized modulated whistlers of the form 

E = ^E(z, <)(x + iy) exp[i(k z — u t)] + c.c, (1) 

where E(z, t) is the slowly varying envelope of the 
whistler electric field, and x and y are the unit vec- 
tors along the x and y axes, respectively, and c.c. stands 
for the complex conjugate. The whistler frequency u 
Ju) ce u) c i) , and the wavenumber kg are related by the 
cold plasma dispersion relation 
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Be 2 : 
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where c is the speed of light in vacuum, w ce = eBo/m e c 
(ujci — eBo/rmc) is the electron (ion) gyrofrequency, 
Upe,o = (47rnoe 2 /TO e ) 1 / 2 is the electron plasma frequency, 
e is the magnitude of the electron charge, m e (mi) is the 
electron mass, and uq is the unperturbed background 
electron number density. 

The dynamics of modulated whistler wavepacket in 
the presence of electron density perturbations associated 
with low-frequency ion-acoustic fluctuations and of the 
nonlinear frequency-shift caused by the magnetic field- 
aligned free streaming of electrons (with the flow speed 
Vez) , is governed by the nonlinear Schrodinger equation 
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u)E = 0, (3) 
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where 



LU = 



^pe + k l c2 



+ k v ez , 



(4) 



and uj 2 e — u> 2 e n e /no is the local plasma frequency in- 
cluding the electron density n e of the plasma slow mo- 
tion. The group velocity and the group dispersion of 
whistlers are 
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Eliminating d z <j> from the governing equations for low- 
frequency density perturbations, we have 



d t v z + \d z v 2 z = - — 
2 mi 



d z \n[ — )- F 

n 



1 



min 



■d z pi. (14) 



The nonlinear Schrodingcr equation for the whistler elec- 
tric field together with the low-frequency equations form 
a closed set for our purposes. 



Dimensionless variables 



and 



1 _ 4^0 ] ( 6 ) 



respectively. 

The equations for the ion motion involved in the low- 
frequency (in comparison with the whistler wave fre- 
quency) ion-acoustic perturbations are 



and 



9 t n 4 + d z (n l v iz ) = 



9tV iz + \d z v 2 iz = -^-d z 

_ f I v4 



m l n i 



(7) 



-d zPi , (8) 



where, for an adiabatic compression in one space dimen- 
sion, the ion pressure is given by p, = p»,o(i»/io) ■ Here, 
the unperturbed ion pressure is denoted by p^o = n T i} 
where T t is the ion temperature. 

The electron dynamics in the plasma slow motion is 
governed by the continuity and momentum equations, 
viz. 



and 



d t n e + d z (n e v ez ) = 



o = ^d z 4>-d z \n(^) F. 



(9) 



(10) 



where T e is the electron temperature, <j> is the ambipolar 
potential, and the low-frequency ponderomotive force of 
electron whistlers is 



In order to investigate numerically the interaction be- 
tween whistlers and large amplitude ion-acoustic per- 
turbations, it is convenient to normalize the governing 
equations into dimensionless units, so that relevant pa- 
rameters can be chosen. We introduce the dimension- 
less variables £ = ujpi.oz / C s , where the sound speed 
is C s = [(T e + 2>Tj)/m iY/\ r = u pi flt, N = n/no, 
u = v z /C s and £ = E / ^ -Aim^T,, + 3Tj); the only free di- 
mensionless parameters of the system are f2 c = Lu ce /u) p i t o, 
k = koc/ajpefi, n = Ti/T e and V g = v g /C s . The normal- 
ized system of equations are of the form 



a. 



-8 6 (Nu), 
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dp 



lnN + l.brjN 2 
1 + 377 



+ a\£\ 2 



(15) 



(16) 



and 



d T £ 



-VgO^ + 1 
1 1 



Pd\,£ 



(17) 



1 + K 2 N + K 2 (1 + K 2 ) 2 \ g 

where the constants are a = (1 + n 2 ) 2 uu 2 efi /uj 2 e K 2 and 
P = (1 + k 2 )(1 - 3n 2 )V 2 /4:K 2 n c . The sign of the co- 
efficient P, multiplying the dispersive term in Eq. (3), 
depends on n: When n < 1/V3, P is positive and for 
k < l/y/3 we see that P is negative. 



III. SMALL-AMPLITUDE SOLITARY WAVES 
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(11) 



The system of equations is closed by means of quasi- 
neutrality 



rii = n e = n, 



(12) 



which is justified if luq < lu cc , is fulfilled with some margin. 
The continuity equations for the electrons and ions give 

Viz = v ez =v z , so that 



d t n + d z (nv z ) = 0. 



(13) 



In the small-amplitude limit, viz. N = 1 + N\, u = u\, 
where N x , u x < 1, Eqs. (l)-(3) yield 



and 



d T Ni — —d^ui, 
dr [ Ul -—\£\ 2 



d T £ = -V g d^£ 



+ i 



d^-m + a^i 2 ), 



Pd^+iN,-^ 



VJ (l + « 2 ) 2 



(18) 
(19) 



(20) 
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where the only nonlinearity kept is the ponderomotive 
force terms involving \£\ 2 . It is important to remember 
that our nonlinear Schrodinger equation for the whistler 
field is based on a Taylor expansion of the dispersion rela- 
tion for whistler waves around a wavenumber ko. Thus, 
this model is only accurate for wave envelopes moving 
with speeds close to the group speed V g , and other speeds 
of the wave envelopes may give unphysical results. Here, 
we look for whistler envelope solitary waves moving with 
the group speed V g , so that Ni and U\ depends only on 
X = £, — Vg T , while the electric field envelope is assumed 
to be of the form £ = W(x) exp(iftr — ik£), where W is a 
real- valued function of one argument. Using the bound- 
ary conditions Ni = 0, ut = and W — at |£| = oo, we 
have k = 0, Ni = -W 2 a/(1 - Vg) and ui = V g Ni. We 
here note that subsonic (V g < 1) solitary waves are char- 
acterized by a density cavity while supersonic (V g > 1) 
envelope solitary waves are characterized by a density 
hump. The system of equations (4) to (6) is then re- 
duced to the cubic Schrodinger equation 




Pd 2 xx W + QW 3 - flW = 0, 



(21) 



where Q = afl c K 2 / (1+k 2 )(1 — V 2 ). Localized solutions of 
Eq. (7) only exist if the product PQ is positive. We note 
that P > (P < 0) when the whistler frequency ojq < 



FIG. 1: Simulations of solitary whistler waves, with its 
associated electric field envelope \£\ (left panels) and den- 
sity perturbation <5iV = N — 1 (right panels). Parame- 
ters are: V g = 640, k = 1.15, fi c = w ce /kV,o = 7.05, 
mi/m e — 73400 (argon ions) and r\ — 0.05. For the initial 
condition, we used the small-amplitude solitary wave solu- 
tion £ = W = (2fi/Q) 1/2 sech[(n/P) 1/2 (C - 2.5 x 10 5 )], with 
/4 (u > Wce/4), and that Q < (Q > 0) when \V g \ < q = _ .025 and P = -7.57 x 10 4 for the given parameters, 

and — —0.5 x 10~ 4 (upper panels) and Q = — 1 x 10~ 4 (lower 
panels). For the density and velocity, we used N = 1 + Ni 



1 (\V g \ > 1), so in the frequency band where lu < w ce /4 
only subsonic solitary waves, characterized by a localized 
density cavity can exist, while in the frequency band luq > 
w ce /4, only supersonic solitary waves characterized by a 
localized density hump exist. Equation (7) has exact 
solitary wave solutions of the form 



and u = Mi, where Ni 



- Vg) and mi = V g N\. 



1/2 



sech 



ft \ 1/2 
i) « 



V 9 r- 



£o) 



(22) 



where V g and ft and the displacement £o are the three free 
parameters for a given set of physical plasma parameters. 
Finally, we recall that the dispersion relation for the elec- 
tron whistlers used here is valid if ujq > ^uj ce uj C i. For sub- 
sonic whistlers having the group speed v g — C s V g (where 



V g < 1), where v g w 2uJo/ko and u>o fcp 



WC 



have cko/ujpe.o = (C s /c)(uj pe $/uj ce )V g /2 > (m e /m i ) 1 / i . 



IV. NUMERICAL RESULTS 

We have investigated the properties of modulated 
whistler wave packets by solving numerically Eqs. (1)— 
(3). We have here chosen parameters from a recent ex- 
periment, where the formation of localized whistler en- 
velopes have been observed 0- In the experiment, one 
has n = 1.2 x 10 12 cm" 3 and B = 100 G, so that 
^ P e,o = 6.7 x 10 10 s _1 and lu C£ = 1.76 x 10 9 s _1 , respec- 
tively. Hence, oj ce /uj pe Q = 0.026. The frequency of the 



whistler wave is ojq — 2n x 160 x 10 s 



1.0 x 10-V 



so that cuo/^ce ~ 0.57 > 0.25. Thus, the whistlers have 



negative group dispersion. From the dispersion relation 
of whistlers, we have k « 1.15, which gives kg ss 257 m _ . 
The latter corresponds to whistlers with a wavelength 
of 2.4 cm. Furthermore, the whistler group velocity is 
v g = 3.36 x 10 6 m/s. The argon ion-electron plasma 
{mi/m e — 73400) had the temperatures of T e = 10 eV 
and Ti = 0.5 eV, giving the sound speed 5.25 x 10 3 m/s, 
and the normalized group velocity V g = v g /C' s = 640. 
In Fig. 1, we have illustrated localized whistler envelope 
solitons, in which the electric field envelope (left panels) 
is accompanied with a density hump (right panels). We 
notice that the density hump is relatively small, due to 
the large (in comparison with the acoustic speed) group 
velocity of the whistler waves. In Fig. 2, we can see the 
development of a large-amplitude whistler pulse, which 
has been launched in a plasma perturbed by ion-acoustic 
waves, with a density modulation of one percent (see 
the caption of Fig. 2). This simulates, to some extent, 
the experiment by Kostrov et o/., where the density and 
magnetic field were perturbed by a low-frequency con- 
ical refraction wave, giving rise to a modulation of the 
electron whistlers. Here, as in the experiment, we ob- 
serve that a modulated electron whistler pulse (middle 
panel of Fig. 2) develops into isolated solitary electron 
whistler waves (lower panel). We note that the wave- 
length of the whistlers is w 2.5 cm, while the typical 
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FIG. 2: The electric field as a function of normalized space 
£ at the times r = (upper panel), r = 2000 (middle panel) 
and r = 4000 (lower panel). The initial condition was a pulse 
on the form £ = 0.07{1 -tanh[5 cos(tt£/3.1 x 10 5 )]}, while the 
density was perturbed as N = 1 + 0.01 cos(207r£/6.2 x 10 5 ). 
Parameters are: V g = 640, k = 1.15, Q c = u C e/t*)pi,Q = 7.05, 
m,i/m e — 73400, rj = 0.05. The (scaled by C a /u pi fi) wave- 
length of the high-frequency wave is Ao = 1.16 x 10 3 , corre- 
sponding to ~ 2.5 cm in dimensional units. The envelope 
±|£| was plotted together with the real part of the wave 
£exp{iK £), where Ko = k C s /u)pi = 2n 2 n c /(l + K 2 ) 2 V g . We 
see the formation of separate wavepackets of the modulated 
high-frequency wave. 



width of a solitary pulse is A£ ss 3 x 10 4 in the scaled 
length units, corresponding to » 64 cm, so that each 
solitary wave train contains 25 wavelengths of the high- 
frequency whistlers. In one experiment, illustrated in the 
lower panel of Fig 4 in Ref. 0] , one finds that the width 
of the solitary whistler pulse in time is 0.2 fj,s, which with 
the group speed v g = 3.36 x 10 6 m/s gives the width 
~ 60 cm in space of the solitary wave packets, in good 
agreement with our numerical results. From the rela- 
tion Ni = —W 2 a/(1 — Vg) valid for solitary whistlers 
in the small- amplitude limit, and with the amplitude of 
W = \£\ approximately 0.3 seen in the lower panel of 
Fig. 2, we can estimate the relative amplitude of the 
density hump associated with the solitary waves to be 
of the order 10 -3 , i.e. much smaller than the modula- 
tion ~ 10 -2 due to the ion-acoustic waves excited in the 
initial condition. 

Next, we study the properties of subsonic whistler 
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FIG. 3: Simulations of solitary whistler waves, with its as- 
sociated electric field envelope \£\ (left panels) and density 
N (right panels). Parameters are: V g = 0.5, k = 0.2, 
tt c = u>ce/upi,o = 0.072, m;/m e = 30000, f] = 0.1. For the 
initial condition, we used the small-amplitude solitary wave 
solution E = W = (2n/Q) 1/2 sech[(fVP) 1/2 (£ ~ 3000)], with 
Q = 5.5 x 10 5 and P — 10.8 for the given parameters, and 
Q. = 1.5 x 10~ 4 (upper panels) and Q, = 3 x 10~ 4 (lower pan- 
els) . For the density and velocity, we used iV = 1 + Ni and 



ui, where Ni 



-W 2 a/{1 



V g ) and mi 



envelope solitary pulses which have the normalized 
group speed V g — 0.5. Here, the restrictive condition 
ck Q /uj pefi = (C s /c)(w pe fi/w ce )Vg/2 > (me/rrii) 1 ^ re- 
quires somewhat higher values of the plasma temperature 
and ujp e fi/uj ce for their existence. With n%i/m e — 30000, 
we have (me/me) 1 / 4 0.1. We take k = cfco/w pei o = 
0.2, C s = 10 5 m/s (corresponding to T e - 1400 'eV) 
7] = 0.1, and uj pefi /uJce = 2400. Thus, fl c = 0.072 and 
u!o/u ce ps 0.039. For these values of the parameters, there 
exist solitary whistler pulse solutions, which we have dis- 
played in Fig. 3. We have used the exact solution in 
the small-amplitude limit as an initial condition for the 
simulation of the full system of equations (l)-(3). The 
bell-shaped whistler electric field envelope is accompa- 
nied with a large-amplitude plasma density cavity. 



V. DISCUSSION 

We have presented theoretical and simulation studies 
of nonlinearly interacting electron whistlers and arbitrary 
large amplitude ion-acoustic perturbations in a magne- 
tized plasma. For this purpose, we have derived a set 
of equations which describe the spatio-temporal evolu- 
tion of a modulated whistler packet in the present of 
slowly varying plasma density perturbations. The pon- 
dcromotive force of the latter, in turn, modifies the lo- 
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cal plasma density in a self-consistent manner. Numer- 
ical solutions of the governing nonlinear equations re- 
veal that subsonic envelope whistler solitons are char- 
acterized by a bell- shaped whistler electric fields that 
are trapped in self-created density cavity. This happens 
when the whistler wave frequency is smaller than uj ce /A, 
where the waves have positive group dispersion. When 
the whistler wave frequency is larger than w ce /4, one en- 
counters negative group dispersive whistlers and the su- 
personic whistler envelope solitons are characterized by 
a bell-shaped whistler electric fields which create a den- 
sity hump. Modulated whistler wavepackets have indeed 
been observed in a laboratory experiment as well as 
near the plasmapause Q and in the auroral zone Q . Our 
results are in excellent agreement with the experimental 



results 2] , while we think that a multi-dimensional study, 
including channelling of whistler waves in density ducts, 
is required to interpret the observations by Cluster and 
Freja satellites. 
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